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ABSTRACT 

Energy spectra of cosmic ray nuclei boron to iron have been measured 
from 2 GeV/nucleon to beyond 100 GeV/nucleon. The data were obtained using 
an ionization calorimeter flown on a balloon from Palestine, Texas. The 
3450kg payload floated at 7g/cm for almost 24 hours. The results are 
in excellent agreement with those of other workers where overlaps exist. 

The spectra are not consistent with single power lav/s, and demonstrate 
the power of using a single technique sensitive over a large dynamic 
range. The data are consistent with the leaky box model of cosmic ray 
propagation. The boron data indicate that the cosmic ray escape length 
decreases with increasing energy as 4+0.1) U p q go GeV/nucleon. 
Secondary nuclei from iron are also consistent with this dependence. 
Predicted changes in the energy dependence of the ratios of primary 
nuclei 0/C and IR0N/C+0 are also observed. 

**Now with University of Siegen, Federal Republic of Germany 
^ Now with Max Planck Institut fur Aeronomie, Federal Republic of Germany 
* NAS/NRC Resident Research Associate, now with Oxford University, England 
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1 . Introduction 

All nuclei of the periodic table of the elements are present in the cosmic 
radiation and the comparison between the relative abundance of cosmic ray 
nuclei and the "the universal abundance" provides information about their 
origin and propagation through interstellar space. Interest in these studies 
has increased since experimental results indicated that the relative abundance 
of cosmic ray nuclei is energy dependent. These results have led to a very 
lively discussion concerning the understanding of the cosmic ray confinement 
mechanisms in the galaxy and interstellar propagation. These observations 
provide clues in the study of the distribution of path lengths through which 
the particles pass in the galactic containment region. The mean interstellar 
matter traversed can be studied by comparing the ratios of the secondary to 
primary components. Observations show that at high energies the particles 
traverse less interstellar matter than at low energies indicating that the 
cosmic ray source composition is less subject to corrections at high energies. 

For this reason much effort has been devoted in the last few years to the develop- 
ment of high energy cosmic ray experiments, and many nuclear species have now 
been ’observed at energies up to 100 GeV/nucleon. These results were obtained 
employing instruments such as ionization spectrometers (Balasubrahmanyan and 
Ormes, 1973, Saito et al., 1974, and Schmidt et al . , 1976), gas Cherenkov, counters 
(Juliusson, J974 and Caldwell, 1977, Lezniak and Webber, 1978) and magnetic 
spectrometers (Orth et al., 1978). Based on these data, a variety of interpre- 
tations were possible, and there were some differences between different tech- 
niques . 

In this paper the results of the chemical composition of cosmic rays as 
a function of energy in the range of a few GeV/nucleon to some hundreds of 
GeV/nucleon for boron through iron are presented. The experiment combined an 
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ionization spectrometer and a gas Cherenkov counter to perform two different 
and independent energy measurements. This combination also provided for the 
first time an experimental cross-check on the response of both types of detectors 
to heavy nuclei. 

In all these balloon-borne experiments the major constraint in extending 
the energy range is the low flux of cosmic ray nuclei at high energies. This 
requires experiments with large area solid angle factors. In order to maximize 
the statistical significance of the results, the experiment was optimized for 
collection area at the expense of both charge and energy resolution. Much of 
this paper is devoted to a discussion of how the resulting problems are handled 
so that the critical reader can evaluate these results. A price was paid in 
terms of the increased background, which had to be evaluated carefully. 

The experiment was successfully flown in October 1976 from Palestine, Texas, 
at a mean altitude of 7 g/cm^ residual atmosphere, providing a total exposure 
of 15.5 m 2 ster hr. The weight of -the scientific payload was 2653.6 kg and the 
total suspended payload was 3450 kg. 

2. Experimental Apparatus 

A detailed description of the experimental apparatus has been published 
(Arens et al . , 1979). Here, a short description of the major components will 
suffice. Figure 1 shows the experimental configuration. This apparatus 
consisted basically of three parts': 

-1) a gas .Cherenkov counter on top of the Instrument with an energy 
threshold of 16.5 GeV/nucleon, 

2) fitfe organic scintillators (Pilot Y) and a solid Cherenkov counter 

(Pilot 425) which measured charge 'and trajectories, and 

3) an ionization spectrometer or calorimeter for energy determination. 

The sensitive area of the apparatus as defined by the upper scintillator 
was 1.2m x 1.2m and the geometric factor was 0.72m^ ster, which reduced to 
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0.28m 2 ster for those particles which penetrated both the gas Cherenkov 
counter and the spectrometer. For the balloon flight, the instrument was 
placed in a styrofoam housing for thermal' insulation. Due to the large 
size of the instrument no pressurized gondola was used. 

The compromise to launch a large payload without a picture device 
(without a pressurized gondola, no gas filled device such as spark chambers 
or multi wire proportional counters could be used) created a background problem. 

Events with a large number of particles, especially due to air showers at 
large zenith angles produced a background which is unusual for experiments 
of this kind. The method of handling this background will be discussed in 
section 5. 

The combination of the gas Cherenkov counter with the ionization spec- 
trometer in this experiment provided a means of cross checking the response 
characteristics of both detector systems to heavy nuclei. The response of 
the ionization spectrometer to the different heavy cosmic ray nuclei could be 
calibrated with the threshold and response of the gas Cherenkov counter; the 
response distribution of the gas Cherenkov counter could be checked for the 
effects of residual scintillation light by using the spectrometer for energy 
cuts. These points are discussed in section 4 and in more detail in a separate paper 
(Balasubrahmanyan et a!., 1979). 

The properties and the limitations of the position sensing scintillator 
technique for determining the trajectories in this high energy cosmic ray 
experiment are discussed in the instrument paper (Arens et al . , 1979). The 
basic idea behind this technique is simple: one locates the position of inci- 

dence by exploiting the fact that the measured signal depends on the location 
in the scintillator relative to the phototubes. Since the pulse height 
measurements are subject to photoelectron fluctuations, this technique is 
most useful for heavy cosmic ray nuclei. The position of an incident particle 
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could be located to within a few centimeters in the different detectors and 
was adequate for our purposes; however, the resulting inaccuracy in the tra- 
jectories did in fact limit the charge resolution. At high energies the charge 
resolution was further hampered by backscatter effects, and multiple tracks 
could not be detected easily, points which will be discussed further. 

3. Charge Determination 

The charge determination was performed with five independent charge 
detectors, i.e. the four scintillators Six, Sly, S3x, S3y and a solid 
Cherenkov counter. Since the Cherenkov counter was viewed by two sets, each 
of 4 phototubes, a total of six charge signals were available. Zenith angle 
corrections and corrections for the spatial non-uniformity of light collection 
were applied to all the charge detectors. A series of scatter plots were then 
constructed from the flight data in which each scintillator was compared with 
the Cherenkov response. The curves through the charge peaks in these two 
dimensional plots, representing the charge lines, were constructed by a best 
fit technique. By using all the detectors, a six-dimensional distribution 
of data and the corresponding charge line was obtained. 


An individual cosmic ray event was attributed to the charge corresponding 
to the nearest point, measured in units of a., on the charge Tine, i.e. by 
minimizing the expression 


i 2 


X 2 = 2 


( X L - X p) 


i=l 


0 ) 


CH 


Here, the x^ are the coordinates of ftfe position on the charge line to be found, 
and the x^ are the charge coordinates of the event considered. The sum extends 
over all the charge coordinates used and each term is weighted by the corres- 
ponding charge signal variance a?. The charge variable was thus effectively 



6 


treated as continuous. The line position was then projected onto the Cherenkov 

axis for an absolute analog charge measure. The Cherenkov axis was chosen, 

since no saturation in the signals up to iron nuclei was found, thus yielding 

a strictly Z ^-proportional charge measure. The above procedure of determini-ng 

the charge avoided any scale adjustment, since the charge Tine took into 

account the individual response of each detector. The absolute value of the 

minimum also provided a measure of the acceptability of an event. Note, 

2 

however, that this variable is not expected to obey a classical x distribution 
because of the non-Gaussian signal fluctuations and a certain degree of correlation 
among the detector signals. 

In the course of the analysis it was observed that the mean pulse height 
for each charge tended to increase systematically with energy as illustrated in 
Figure 3. The rate of signal increase is more pronounced the closer theddetector 
is located to the spectrometer. This behavior is 

attributed to relatively low energy particles which flow back from the iron 
spectrometer. These particles then stop preferentially in the detector close to 
the spectrometer. ‘The increases in SIX and S1Y, the detectors which are further 
away, are quantitatively consistent with being caused by 6-ray production between 
1 and 20 Gev/nucleon (Yodh, 1977). Because of the backscatter, only SIX, SI Y 
and the two sets of PMTs from the solid Cherenkov-detector have been used, and 
the x^ values should correspond to 3 degrees of freedom. Fig. 2 shows a charge 
histogram for particles in the energy range 20 to 40 GeV/nucleon obtained by 
accepting all events with x" - values smaller than 10 and a comparable background 
histogram for x 2 values between 36 and 100. (A discussion of handling events 
in the x 2 range 10-35 is deferred until section 5.) The more plentiful species 
(B, C, 0, Ne, Mg, Si, Fe) stand out conspicuously and the 
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charge resolution is about 0.4 charge units around oxygen and 1.6 charge units 
around, iron. The detector resolution (a) was about 10% independent of charge and 
was limited by the accuracy with which the trajectories of an incident particle 
could be determined., Because individual element resolution has not been obtained, 
the elements have been grouped together as indicated in the figure. Because 
nitrogen lies between the abundant carbon and oxygen and does not clearly stand out, 
spectra for this nucleus should be considered cautiously. The boron peak appears 
to stand out in all the plots and so those results are included. The separation 
of boron from the abundant background will be discussed later. Note how steeply 
the background decreases as a function of increasing charge. 

The' signal resolutions a. ‘s in equation 1 taken for the four detectors SIX, 

SI Y and the two sets of Cherenkov tubes were respectively 7, 5, 7 and 9 percent. 

They were obtained by a subset of the carbon and oxygen nuclei which penetrated 
the central regions of the detectors. We found them to be independent of charge. 

We assumed that they are independent of energy. This might lead to a slight 
underestimation of a*s at high energies, but as will be shown later, an indepen- 
dent decision is made about the background correction for each charge and energy 
bin, taking account of any variation. An incident iron nucleus which strips off 
an alpha particle at the bottom of SI Y and hence registers as a chromium in the 
Cherenkov detector would, in the absence of other fluctuations give a x 2 value 
of 9 and lies at the limit of our detection. An iron that fragments into a calcium 
on the other hand, would give a x 2 value greater than 50 and is easily rejected. 

An allowance is made for the undetected interactions when the correction is made 
for interactions in the telescope. 

4. Energy Determination 

The principal detector for determining the energy of the cosmic ray nuclei 
was the ionization spectrometer (Fig. 1). A primary particle converts a sub- 
stantial part of its initial kinetic energy into a nuclear electromagnetic cascade. 
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The cascade builds up by a series of nuclear interactions in the spectrometer. 

Energy is measured by sampling the resultant electromagnetic cascades at 

various depths in the absorber. In this experiment the iron slab scintillator 

sandwich was divided into three modules. • Each module was viewed from two 

opposing sides by phototubes at the ends of air light guides with white walls. 

2 

A vertically incident cosmic ray encounters a total of 170 g/cm iron plus 
5.9 g/cm^ scintillator material, which corresponds to 1.44 proton and about 6 
iron nuclear interaction lengths. 

For an infinitely deep calorimeter the SUM signal from all the scintillator 
samples is proportional to the kinetic energy of the incident nuclei no matter 
where the first interaction within the calorimeter takes place. All the energy 
is absorbed and the signal distribution would be Gaussian in shape and have a 

resolution of about 10%. However, in this finite thickness device, the fraction of 

\ 

the primary energy deposited is dependent upon the location and characteristics 
of the first interaction, resulting in a broad SUM signal distribution. A 
thorough discussion concerning this effect on the shape of a measured cosmic ray 
spectrum is given in an earlier paper, Jones et al , (1977). In this analysis 
we utilize the fact that the cascade develops differently depending on the location 
of the first interaction in order to improve the resolution 

The gas Cherenkov counter has been used to provide a calibration for this 
purpose. Fig. 4 shows the response of the gas Cherenkov counter for oxygen nuclei 
as a function of the calorimeter signal. No restrictions have been placed on the 
calorimeter signals but data have been selected to have trajectories passing: 
through the gas Cherenkov counter and to have consistent signals in the two banks 
of PMT's. At low calorimeter signal levels the .finite response in the gas Cherenkov 
counter is due to the scintillation in the Freon gas. As the energy of particles 
can be identified in this experiment using the calorimeter, this scintillation 
level is determined explicitly. Consequently the threshold of the gas Cherenkov 
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counter can be determined quite accurately minimizing the uncertainty due to the 
scintillation contribution of the gas Cherenkov counter. The curve, shown in 
Fig. 4 is based upon the energy determined from Equation 2, described below.. 
Similar plots are obtained for other nuclei. From these plots the charge- 
independence of the calorimeter signal (on a per nucleon basis) corresponding, 
to the threshold of the gas Cherenkov counter is demonstrated as seen in the- 

insert in Fig. 4. 

Plots of the gas Cherenkov response as a function of the calorimeter signal 
can also be made for given Fel signals. In this manner, the correlation, of 
the calorimeter SUM signal with the Fel signal can be found. This correlation is 
indicated by the data points in Fig. 5. In addition to this experimental check, 
Monte Carlo simulations (Jones, 1976 and Jones et al , 1977) have been used, to: 
deduce the correlation for other energies. An empirical fit was found which- is 
wall represented by the following expression: 


(Fe l) 1 


(GeV/nucleon) . 


The lines in Fig. 5 are for constant energy as computed by this formula... The 
term SUM 1 ' 36 has the effect of allowing for the increasing fraction of energy 
escaping out the- bottom of the calorimeter as the energy increases, and the term- 
Fei -0.38 compensates for what happens in the first interaction. Note that when 
the SUM and Fel signals are expressed in particles per incident nucleon,. this 


expression is independent of charge. The Monte Carlo calculations and the gas 
Cherenkov data gave consistent values for the constant which represents the 
energy calibration. Thus, at 16.5 GeV/nucleon our energy scale has been deter- 
mined- to within 5% and is independent of charge. At other energies, since we 


depend on the Monte Carlo simulation, an increased energy scale error is allowed" 
because of the distance from the gas Cherenkov calibration point'. These esti- 


mated errors are summarized in Fig. 6. 
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An example of how equation 2 fits the Monte Carlo simulation is shown 
in Fig. 7 for 100 GeV/nucleon carbon. The clump of particles in the lower 
left hand corner of the plot have penetrated deeply into the calorimeter’ and 
are lost for energy determination. A correction for this effect is discussed 
below. The simulations show some Z dependent effects at low energies (1-3- GeV/ 
nucleon) which may indicate slight problems with the Monte Carlo program at 
energies where ionization loss effects are important. Therefore only data' 
for energies E >_ 2.5 GeV/nucleon are analyzed. 

Observed energy distributions for Monte Carlo simulated data are shown 
in Fig. 8. Simulated monoenergetic incident carbon nuclei of 30 GeV/nucleon 
are analyzed by two different methods. The solid histogram shows the distri- 
bution obtained by taking the SUM signal and the dashed histogram shows the 
improvement obtained by applying the equation (2). The dashed histogram peaks 
more sharply around 30 GeV/nucleon although both distributions extend to near 
zero for those particles which penetrated deeply before interacting. 

The Monte Carlo events provide a series of distributions which can then be 
used to deconvolve the energy spectra. Because of the. power law nature of the 
energy -spectra the slopes are not very sensitive to these deconvolutions. This 
is true as long as the shape of the response distribution varies slowly with 
energy (Jones et al., 1977). On the other hand, the intensities must be corrected 
for the spillage from bin to bin. For example, particles which penetrate deeply 
into the calorimeter with little or no interaction are lost. They appear as low 
energy particles, and because of the steeply falling cosmic ray spectrum they 
have a statistically significant effect only when they fall in the next lower 
energy bin. Particles which fall into the next higher energy bin are few in 
number but their effect is enhanced by the cosmic ray spectrum. In order to 
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deduce the correct fluxes Monte Carlo simulated distributions for carbon, 
silicon and iron nuclei at 3, 6, 17, 30 and 100 GeV/nucleon were used. Energy 
bins were overlaid on these plots, and the events lost and those going into 
adjacent bins were counted. By allowing both for the events lost and gained, 
appropriately weighted by the spectrum, a correction was derived.- Fig. 9 shows 
the multiplicative correction factor (CE) derived. At low energies nuclei stop 
by ionization loss so the correction factor is 1, but at higher energies particles 
can penetrate deeply without interacting or by suffering low multiplicity 
interactions. The net result is a slight flattening of the observed spectra. 
Uncertainties in these deconvolution corrections are based on the statistical 
significance of the Monte Carlo runs and are indicated in Fig. 9. This error has 
been included in the quoted flux and dominates at all but the highest energies 
where statistical errors dominate. In taking ratios of intensities the effect of 
these convolution corrections tend to cancel. 

5. Background 

Because of the large area of the experiment and the lack of a trajectory 
picture device, a substantial background of multiple particles exists, -apparently 
due to air showers, mostly from large zenith angles. Most of these 'events lie 
at large x 2 values (see equation 1 and the discussion of the charge analysis) 
indicating that thejr pulse height were inconsistent in the SIX, S1Y and the ■ 
Cherenkov detectors. However, the low x 2 tail of the distribution can simulate 
good events. The background falls off rapidly with increasing charge but because 
it has an "apparent" spectrum which is considerably flatter than that of "real" 
•events, is troublesome for low Z particles at high energies. A background event 
usually gives: a trajectory (which measures the "center of gravity" of particles) 
which is reasonable and so this cannot be used to reject these events. In summary. 
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the relative amount of background is both charge and energy dependent. The 
method of determining the flux of particles in a given charge and energy bin 
is illustrated in Figure 10. The open circles indicate the x 2 distribution 
for all carbon nuclei in the bin centered at 15.8 GeV/nucleon. Erelow x 2 = 

10, the distribution is dominated by good events with resolved charges. At'high 
X 2 , beyond 36, the distributions show no evidence of charge resolution and 
the distribution is clearly dominated by background. The problem is to determine 
what to do in the transition between where there is an admixture of background 
and good events, and to find an effective cut that properly: estimates the 
background. To do this, a x 2 distribution consisting predominantly of background 
was constructed, based upon the shape of the cascades in the calorimeter. 

For legitimate downward moving high energy particles, cascades should build 
up with increasing depth (e.g. Fe 3 signal > Fe 2 > Fe 1). It was observed 
that events for which this condition was violated had indeed mainly high x 2 
values. Therefore, by selecting events for which a major fraction of the energy 
was deposited in Fe 1, ax 2 distribution was formed for background events for 
each charge and energy bin. The example shown in Figure 10 by the solid circles 
is typical, and the shapes (for background) were nearly independent of energy. 

The background distribution is then normalized to the all particle distribution 
at large x 2 values and the flux of good events is found by a subtraction of two 
distributions. An effective cut can also be defined which gives the same 
intensity. At lower energies, the peak to valley- ratios are pronounced (10:1), 
and the cuts are obvious and typically occur at x 2 values of 30. A sequence 
of examples is shown in Figure 11 for boron, carbon and iron secondary nuclei 
in the energy intervals indicated. The error in this procedure is comparable to 
the statistical error and has been included in the quoted errors. 
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6. Interaction Corrections 

2 

A correction is made for the interactions in the 2.57 g/cm of As, and the 
o 

1.28 g/cm of scintillators at the top of the instrument. The interaction -mean 
free paths and the corrections used are given in Table 1. The correction factor 
is adjusted for the fraction of interactions which can be detected by the x 2 
■techniques used. For example, in the case of iron secondaries (17 £ Z £ 20) the 
correction factor would be 1.43, but since only 2/3 of the interactions are detected, 
this reduces to (1 + 2/3 x .43) = 1.29. The factor 2/3 -comes from the partial cross 
sections for interactions which go undetected, i.e. 1 in 3 interacting nuclei do not 
chanqe charge by enough to raise their x 2 sufficiently to be rejected. 

To correct the fluxes obtained to the top of the atmosphere, a full propa- 
gation analysis has been done and growth curves derived. The mean free paths 
used are also shown in Table 1. For boron .and iron secondaries, a range of .values 
is given, because the correction depends on the abundance ratio (to carbon and iron 
respectively) observed; the smaller the observed abundances, the larger the fraction 
which has been produced in the atmosphere. The quoted errors indicate the range of 
■values used for these data. This correction becomes very important at high energies 
where the interstellar material traversed is less (at 100 GeV/nuc = 1 g/cm ) and 
more than half of the observed flux is atmospherically produced. 

7. Results and Conclusions 

The derived fluxes of all the compqnents at the top of the atmosphere are 

2 5 

given in table 2. Results for oxygen and iron nuclei, multiplied by E ’ power 
to emphasize the spectral structure, are shown in Figure 12 along with some of 
the recent data from other groups. The agreement between the different measure- 
ments is quite good and generally within the quoted errors. The oxygen data of 
Lezniak and Webber (1978) seem to be slightly lower than those of the other groups, 
e.g. the Orth et al . (1978) data and those of the Chicago group (Juliusson, 1974 
and Caldwell , 1977) . 
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The spectra derived from this experiment are shown in figure 13, again 
multiplied by E^' power. The spectra are given as functions of kinetic energy. 

The oxygen data span almost 3 decades in energy (corresponding to 7 decades 
in intensity). The low energy data, (2-20 GeV/nucleon) when plotted as a function 
of total energy x Ey* , are quite flat; but the data at high energy indicates 
that there is a gradual steepening of the spectra. For example, between 2 and 

o p 

20 GeV/nucleon the iron data is fitted by E^ ’ while the highest energy data 
is consistent with E^, * . The results agree with spectra derived using a • 
thicker calorimeter (Balasubrahmanyan and Ormes, 1973) and with gas Cherenkov 
counters (Juliusson, 1974, Caldwell, 1977) when allowance is made for the different 
energy coverage of the different experiments. It should be noted that the 
statistics are not sufficient to rule out the iron spectrum Goodman et al. (1979) 
require to explain their observations of delayed particles at mountain altitudes. 
Their result at 2-20 TeV/nucleon indicate an iron spectrum as flat as E“2.36+.06 
from 10 GeV/nucleon onwards. 

In Figure 14 the experimental B/C+0 ratio is shown along with the best fit 
power law dependence of the mean escape length, = 9 g/cm2 (E/GeV/nuc)"^ 0 * 4 — 

The calculation is based upon an exponential distribution of vacuum path lengths 
in the phenomenological leaky box model for cosmic ray propagation. The inter- 
action mean free paths used for the propagation calculation, including the 
effect of 10% He, are given in the last column of Table 1. (See Ormes and Freier 
1978 and references therein for a discussion of the propagation model). This 
mean escape length has been allowed to vary with energy in the manner indicated 
in order to fit to the data in the energy range of this experiment. This energy 
dependence of the mean escape length is consistent within errors of those previously 
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derived (Caldwell, 1977, Lezrriak and Webber 1978 and Orth et al. 1978) Because 
of the competition between escape and interaction loss mechanisms, ratios between 
primary nuclei also vary in this model. In Figure 15 and 16 the 0/C and the iron/ 
C+0 ratios are shown, again compared with the expected ratios as a function of 
energy. The X has the same form as fitted to the B/C+0 data, and the source 
ratios are taken from Shapiro et al. (1975). A 10% higher source abundance of 0 
would fit the 0/C ratio and the IRON/C+O ratios better. The larger source ratio 
of iron derived by Lezniak and Webber (1978) is not indicated by these data. 

Finally in Figure 17 the ratio of iron secondaries (17<Z<25) to iron nuclei 
is shown along with the results from other workers. Our high energy data tend to 
fall below the predicted ratios, especially in the energy range 3-30 GeV/nucleon. 

The predicted ratios have been calculated based upon a cross section for frag- 
mentation of iron into (21<Z<25) of 411 mb. The higher ratio observed around 
1 GeV/nuc by Lezniak and Webber (1978) probably reflects the increase in the cross • 
section towards lower energy. The sub-iron nuclei can be used to probe the shape 
of the cosmic ray path length distribution. Various workers (Lezniak & Webber, 

1979, Garcia Munoz et al., 1977) studying cosmic rays below 1 GeV/nucleon have 
suggested that the sub-iron abundances are large enough that the path length dis- 
tribution must be truncated below 1 g/cm 2 (Shapiro et al., 1973). Since the cross 
sections and especially their energy dependence in the range 1 to 10 GeV/nucleon 
are not adequately known, it is difficult to draw any conclusions about the cosmic 
ray path length distribution based upon these data. However, since truncation of 
the cosmic ray path length distribution will increase the number of iron secondaries 
predicted, it is difficult to see how our data can be consistent with any trunca- 
tion. Further measurements of the cosmic ray intensity above the atmosphere and 
better determination of the energy dependence of the cross sections of the energy 
range 1-10 GeV/nucleon will be required to resolve this question. 
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In conclusion, the data are all consistent with a quite simple propagation 
model in which the escape of cosmic rays from the galaxy is energy dependent. At 
high energies (>50 GeV/nuc) the uncertainties in the atmospheric corrections are 
important and indicate the need for satellite measurements. At around 100 GeV/ 
nucleon the cosmic rays observed at earth are almost an order of magnitude younger 
than those at 1 GeV/nucleon and more nearly reflect the source composition. When 
high statistics observations of the nuclei in the charge range 9<Z<25 are made 
above the atmosphere at 100 GeV/nucleon the source composition of these rare 
components can be determined. 
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TABLE 1 


CHARGE 

MFP* 

IN 

ALUMINUM 

MFP* 

in 

SCINTILLATOR 

Correction for 
interacting 
events rejected 

MFP* 

in 

AIR 

Atmospheric* 

Correction 

MFP* 

in 

ISM' 

(10% He) 

5 

43.2 

21.3 

1.20 

27.5 

,76-. 96. 

6.6 

6 

40.1 

20.0 

1.22 

26.2 

1.25 

6.2 

7 

37.2 

18. 6 

1.23 

24.0 

1.22 

5,6 

8 

34.7 

17.5 

1.26 

22.7 

1.32 

5.1 

10 

30.8 

15.7 

1.29 

20.8 

1.23 

4.4 

12 

27.8 

14.4 

1.33, 

10.4 

1.34 

3.9 

14 

25.6 

13,3 

1.36 

18.2 

1.40 

3,52 

1 0 <Z<1 6 



1.33 


1.34 


1 7<Z<_20 



1,29 


0,86-1.00 


21<Z<25 



1.31 


0.62-0.93 


26 

19,8 

9.3 

1,43 

13,5 

1.44 

2.2 


2 

*MFP = mean free path for interaction in g/em - 

+ p 

At mean slant depth of 7.5 g/cm 



TABLE 2 


Energy 

GeV/Nucleon 

Boron 

Carbon 

2.5 

0.17+o03 

0 . 68 j> 1 

4.0 

(7.7+1 .2) -2 

0 . 27+ . 04 

6 . 3 

(2.0+.3J-2 

(9.4+1 .4)-2 

lOo 

(5.6+.D-3 

(3.5+o5)-2 

16. 

(1.8+o3)-3 

( 1 <, 2 + o 2)-2 

25. 

(7. 9+1 .2 )-4 

(3,,7+.5)-3 

40 o 

(1 .5+. 5 )-4 

(l.l+.2)-3 

63. 

(4o0j4 .4 )-5 

(2. 5+.5)~4 

100 . 

(l + «4 )-5 

(8.0t.D-5 

160. 


(2„9+o4)-5 

250. 


(4.4+1 .3)-6 

400. 


( 1 » 8 +o 8)-6 

630. 


(6ol+o4)-7 


p 

(Fluxes in Particles/m sr sec GeV/nuc) 


(X+AX)-N 

= (X+AX) x 

10 " N 

Nitrogen 

Oxygen 

1 0<Z<1 6 

0. 25+« 04 

0.66+.09 

0o48+.07 

0 . 11+.02 

Oo 29+o 04 

0.22+o03 

(3.4+.5)-2 

0.10+.013 

(7.9+1.2)-2 

( 1 . 2 +o 2)-2 

(3.8+ 0 5)-2 

(2.8+.4)-2 

(3.7+.6)-3 

(1.3+o2)-2 

(.99+. 15) -2 

(Ul + o2)-3 

(4o5+o7)-3 

(2.9+.4)-3 

(2.7+.5)-4 

(1.3+.2)-3 

(o95+.15)-3 

(7.4+2. )-5 

(3o3+o5)-4 

(3.0+.5)-4 

(1 0 8+o6)-5 

(1o3+o3)-4 

( ,91+.2)-4 

(3.9+2 0 )-6 

(lo8+.4)-5 

(2.4+o5)-5 

(4 0 + 4)-7 

(5.9+1 o3)-6 

(. 8 I+ 0 2) -5 


(3o2+lo0)-6 

(1.7+;7)-6 


(1o2+o5)-6 

(5.4+3)-7 


(2.2+1.3)~7 

(lo2+1.2)-7 


1 7<Z<20 

2 1 <Z<_25 

26<Z<30 

(2.6+.4)-2 

(1 .5+.- 3) -2 

(5.8+0.9)-2 

(1.2+.2)-2 

(8.2+1 .3)-3 

(3.0+o4)-2 

(4.5+.7)-3 

(3.0+.5)-3 

(1.3+.2)-2 

(l.,7+.3)-3 

(1.0+.2)-3 

(4.2+.7)-3 

(5.0+.8)-4 

(3.5+.6)-4 

(1.6+.3)-3 

(1 .4+. 3) -4 

(1 .l + .2)-4 

(4.8+.9)-4 

(5.3+1 .1 )-5 

(3.7+.8)-5 

(1.9+.4)-4 

(1.2+.4)-5 

(1.0+.3)-5 

(3.9+1 .2) -5 

(6. 5+2. 3) -6 

(2.7+1 .1 )-6 

(4.6+3. )-6 

(1.3+.8>6 

(7,6+4)-7 

(6. 9+3) -6 

(2o6+2)-6 


1000 » 



FIGURE CAPTIONS' 


Schematic diagram of the instrument. 

Charge histograms for particles between 20 and 40 GeV/nucleon. The 

upper curve is for good events, and the lower curve is for background 

events as selected by the consistency of measurements in the charge detectors. 

Mean response of charge detectors as a function of energy. S^ and S~ are 

ox dy 

detecting backscattered energy from the calorimeter and hence are not used 
in the charge analysis. 

A cross plot of the gas Cherenkov detector and the calorimeter for oxygen 
nuclei. Such plots are used to check the calibration of the calorimeter in 
the Monte Carlo simulation. The insert shows this calibration for different 
incident nuclei, confirming its independence of charge. 

This curve shows lines of constant energy on a SUM vs. Fe 1 crossplot. 

The background region is where events fall which have cascades which develop 
too fast in the calorimeter. The data points indicate the calibration from 
the gas Cherenkov threshold based on in-flight oxygen nuclei. 

The estimated energy uncertainty as a function of energy. 

Calorimeter response cross plot from the Monte Carlo simulations of 
100 GeV/nucleon carbon nuclei. The solid curve represents the calibration 
curve for 100 GeV/nucleon particles. 

The solid curve shows the calorimeter SUM response distribution for 30 GeV/nucleon 
Monte Carlo simulated carbon events. The dashed curve shows the distribution of 
observed energies when corrected for the Fe 1 pulse height and the energy 
escaping out the bottom using equation 2. 

Multiplication factors which account for energy bin spillover effects. 

Primarily these particles have penetrated too deeply in the device to have 
their observed energy within the resolution bin. 



22 


10. An example (15.8 GeV/nucleon carbon nuclei) to illustrate the procedure 
for subtracting background. The background (filled circles) distribution 
comes from events in the cross hatched region of Figure 5. 

11. Several examples for different charges and energy bins of the x 2 distri- 
butions and the renormalized background distributions. Fluxes are obtained by 
subtracting these two distributions. 

12. The oxygen and iron spectra as dete’rmined from this experiment (solid circles) 

compared with measurements of other values. The spectra have been multi- 
plied by (kinetic energy) to emphasize differences. The intensities 

are in parti cl es/m^ sr GeV/nuc and kinetic energy is in GeV/nuc. The open 
squares are due to Orth et al. (1978), the open circles to Caldwell (1977), 
the crosses to Juliusson (1974), and the vertical lines to Lezniak and 
Webber (1978). 

13. The spectra of various components as determined from this experiment multi- 
plied by e 2*5 (kinetic energy) to emphasize spectral features. The inten- 
sities are in particles/m2 sr GeV/nuc and kinetic energy is in GeV/nuc. 

14. The boron to carbon plus oxygen ratio from these measurements compared to 
mean escape lengths calculated as a function of the energy dependence 
i-ndi cated. 

15. The same calculation (see Figure 14) compared to the 0/C ratio data. The 
source ratios are taken from Shapiro et al., 1975. 

16. The same calculation (see Figure 14) applied to the IRON GR0UP/C+0 ratio data. 

17. The ratio of iron secondary nuclei to iron nuclei compared to the data of 

other workers. The symbols are as follows: G, Garcia-Munoz et al., 1977; 

L, Lezniak and Webber, 1978; C, Caldwell, 1977; J» Juliusson, 1974; 0, Orth 
et al., 1978; and Solid Symbols represent this work. 
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